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FOREWORD 


Series  expansion  formulae  have  been  applied  to  the  Gaussian  form  of  the 
Lagrange  planetary  equations  in  order  to  develop  a  compact  formalism  that  repre¬ 
sents  satellite  drag  perturbations  produced  by  an  exponentially  varying  atmospheric 
density  distribution.  The  compact  result  lends  itself  to  facile  manipulation, 
either  manually  by  the  user  or  automatically  by  the  electronic  computer.  This 
document  has  been  reviewed  and  approved  by  R.  W.  Hill. 
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INTRODUCTION 

In  recent  years,  much  attention  has  been  devoted  to  the  analytic  representa¬ 
tion  of  artificial  satellite  aerodynamic  disturbances  produced  by  the  earth's 
atmosphere*  ^ .  Such  representations  tend  to  generally  address  satellite  orbits 
in  the  small  eccentricity  regime  (e.g.  e  <,-2),  so  that  truncated  expansions 
about  the  eccentricity  can  be  used  to  simplify  the  complicated  perturbation 
equations.  Quite  satisfactory  results  are  produced  in  this  way  for  a  variety  of 
assumed  atmospheric  density  structures. 

This  report  has  been  prepared  to  demonstrate  that  when  one  is  concerned  with 
an  exponentially  varying  atmospheric  density  distribution,  the  aerodynamic  pertur¬ 
bation  equations  can  be  written  in  a  compact  form  without  resorting  to  series 
truncation.  This  formalism  lends  itself  to  computational  expediency,  especially 
for  applications  using  the  electronic  computer,  and  allows  the  user  to  readily 
inflate  the  expressions  to  any  desired  order  in  the  orbital  eccentricity. 

The  following  sections  of  this  report  discuss  in  detail  the  derivation  of 
this  compact  formalism.  The  Gaussian  form  of  the  planetary  equations  and  their 
preliminary  manipulation  into  a  more  utilizable  configuration  are  presented  in 
the  next  section.  The  operations  used  in  the  expansion  process  of  the  derivation 
are  then  discussed  and,  in  the  following  section,  are  applied  to  the  reconfigured 
planetary  equations  to  yield  the  desired  results. 
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THE  VARIATION  OF  THE  ELEMENTS 


The  Gaussian  form  of  the  equations  describing  the  variation  of  the  Keplerian 
element  set  a,  e,  i,  to,  and  ft  are  well  known  and  given  by^ 


a  =  2a  (  — 


(?h 


•  /Vl  -  sin  E\  ,,  .21  cos  E 

e  =  - - - j  Fr+-It—-—e 


(l)-  =  (JLC2S_u)  Fc 


ft  =  r  \ln-U  )  F 

h  \sin  i  C 


/cos  E  +  e\  ,  Vl  ~  e"  /  sin  E  \  „  A 
"  -  ^ - j  Fr  +  2  - ll  -  e~ cos  E>  FI  " 


ft  cos  l 


where  u  is  the  true  argument  of  latitude  defined  as 

u  =  e  +  to  (6) 

0  is  the  true  anomaly,  and  E  is  the  eccentric  anomaly.  The  quantities  r,  v,  and 
h  are  the  satellite  radial  distance,  speed,  and  angular  momentum,  respectively, 
and  are  given  by  the  following  expressions: 


-  e 

cos 

E) 

u 

1  (L 

4-  e 

cos 
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1  ll 

-  e 

cos 
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and 


h  = 


[via(l  -  e2)] 


h 


(9) 


where  p  is  the  earth's  gravitational  constant. 

The  acceleration  components  F^(i  =  R,  I,  C)  appearing  in  the  Gaussian 
equations  are  defined  by  the  vector  equation 


_  A  A  A 
F  =  FRR  +  FXI  +  FCC 


(103 


where 


A 

C  = 


r  x  r 


r  x  r 


A 

I 


and 


(11] 


(12) 


A 

R  = 


A  A 
I  x  C 


(13) 


Here  r*  and  r  are  the  satellite  position  and  velocity  vectors. 

The  drag  force  Fp  per  unit  satellite  mass  can  be  well  represented  by  the 
equation 


SA  - 

~2m  ’Va 


(i4: 


where  CD  is  the  drag  coefficient,  A  is  the  effective  cross  sectional  area  of  the 
satellite  normal  to  the  direction  of  motion,  m  is  the  satellite  mass,  p  is  the 
local  atmospheric  density,  and  V  is  the  relative  velocity  given  by 
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this  expression,  V  is  the  satellite  velocity  and  V  is  the  atmospheric  rotationaJ 

a 


locity.  To  a  very  good  approximation  the  atmospheric  drag  acceleration  components 
n  be  expressed  as7 

1-  2 

C_A  /  5\2  ,,  „v2 

„  D  a  „  .  .  _  (1  -  e  cos  E)  •  . , , . 

F„  =  -  a  ( —  I  epVw  cos  \  sin  E  - ^  E  (16) 


2m 


(1  +  e  cos  E) 


CrjA  ,  i.  . 

FT  =  -  — —  aapV(l  +  e  cos  E)2(l  -  e  cos  E) 2  E 
I  z  m 


i. 


CdA  /5\‘  2 

=  -  —  —  a  I  —  I  pVme  sin  i(l  -  e  cos  E)  cos  u  E 


iere  is  the  earth's  angular  rotation  rati  and 


1  - 


3, ,  2.-1  2 

a  (1  ~  e  ) 


We  cos  t 


(17) 


(18) 


(19) 


Assume  now  that  the  atmospheric  density  varies  exponentially  with  r  and  has 
le  form 


P  =  Ppexp[6(rp  -  r)] 


P 


ppexp£-8ae  + 


Bae  cos  E] 


(20) 


iere  t  Is  the  local  atmospheric  density  at  perigee  and  8  is  the  inverse  density 
P 

•ale  height  that  is  assumed  to  be  constant.  Using  Equations  (7),  (8),  (9),  (16), 
.7),  (18),  and  (20)  along  with  the  identities 


2 

cos*  u 


sin  w  +  cos  2ii 


cos  E  -  e 
1  -  e  cos 


sin 


sin  E(cos  E  -  e) 
(1  -  e  cos  E)7 


(21) 
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sin 


in  Equations 
the  (.a  iss  i  at 
following  ft' 

da 

dE 


de 

dE 


ii 

dE 


u  cos  u  =  -sin  u,  cos  to  +  sin  2to 


cos  E  -  e 


+  i  1  e“  )  COS  2(.> 


1  -  e  cos  E 
sin  E  (cos  E  -  e) 


(1  -  c  cos  E)‘ 


1  1)  tmuu^h  (5),  one  finds  that  after  much  algebraic  manipulation 
equations  tor  the  variation  of  the  elements  can  be  cast  into  the 

riris  : 


- •a“exp[bae  cos  E](l  -  e^  cos^  E)  (1  +  2e  cos  E  +  e^  cos^  E)  (: 
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2  2 
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+  2a  [cos  E  +  e  cos^  E^j 
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2  2 

exp[8ae  cos  E^(l  -  e  cos  E) 


!2  2  2 
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A  2  3  2  3 
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2  2  2  A  2 

-  e  (cos  2w  +  e  sin  to)  cos  E  +  (1  -  e  )  sin  2w  sin  E 
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d n 
dE 


2  we 


u(l  -  e  ). 


2  2 

exp[f3ae  cos  E](l  -  e  cos  E)  • 


jsin  w  cos  u)[(2e^  -  1)  -  2e  cos  E  +  2(1  -  e^) 


4,  2  p 

cos  E 


+  2e^  cos^  E  +  e^(e^  -  2)  cos4  e]  -  (1  -  e^)  cos  2to  sin  E  • 


T  p  2  2  „  ,  2  3 

[e  -  cos  E  -  e  cos  E  +  e  cos  EJ> 


(26) 


_c 


dco  ,  d£i  •  Aa  p  irT  / 1  2  2  .  r 

—  +  cos  i  =  —  exp  |_6ae  cos  E J  ( 1  -  e  cos  E)  sin  E 
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2  2  2 
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(1  +  e  cos  E) 
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re 


A  = 


Cr  A 
D 


m  o 


ap  exp  ^-6ae] 


(28) 


Finally,  in  order  to  obtain  the  changes  produced  in  the  orbital  elements 
atmospheric  drag  during  one  orbital  revolution.  Equations  (23)  through  (27) 
be  integrated  with  respect  to  E  over  2n  radians;  i.e.. 


(29) 
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INTEGRAND  EXPANSION  AND  THE  J  ,  FUNCTION 

n,k 

In  this  section,  a  novel  method  will  he  used  to  provide  expansion  formulae 
for  Equation  (29).  As  can  be  seen,  the  right-hand  side  of  this  equation  is 
composed  of  sums  of  integrals  of  the  form 

211 

/2  2  5 

exp [ Lae  cos  F.](l  -  e  cosz  E)  cos’1  E  dE  (3C 

u 

The  radical  appearing  in  t.he  integrand  of  can  readily  be  expanded  using 
the  Legendre  polynomial  generating  function^ 


(1  -  2tx  +  t  ) 


c\<x) 


if  one  lets 


x  =  0 


t  =  i  e  c 


•os  E  (i  =  -yTT) 


The  P^(x)  In  Equation  (31)  is  the  Legendre  polynomial 
When  Equations  (32)  and  (33)  are  used,  then 


n  ,  2  2 

cos  L(i  -  e  cos  E) 


?  f 


2  k+n 

cos  t 
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in  Equation  (34),  one  finds  that  Equation  (30)  may  be  written  as 


R  (Bae)  =  /  e  J  (Bae) 

n  /j  n  ,k 


In  this  expression. 


,  rsaex  =  +_n)!  ,  /2k 

n,k(Sae)  24k  +  n  -  I 


2 

VeS 


(k  -  j  + 


I 2 j (Bae)  +  2A 


2  r 

n,0  E  f 
j=0  L 


(k  +  j  +  ILyJ-)  !  (k  -  j  +  -1) 


where 


I2j+l(6ae) 


1  for  even  n 


0  for  odd  n 


1  for  odd  n 


0  for  even  n 


and  1  (j'ae)  is  a  Bessel  function  of  the  fir  s'- 1  kind  and  inncinarv  uruurent 
m 


defined  by 


l-ae)  =  f  exp  f.-ae  cos  !■]  cos  (rl')d 


C|  eg 


NSWC  TR  84-37 


THE  CONSUMPTIVE  DRAG  PERTURBATION  EQUATIONS 


It  should  be  noted  that  those  portions  of  the  right-hand  sides  of 
Equations  (25),  (26),  and  (27)  that  contain  sin  E  multipliers,  will  vanish 
from  the  right-hand  side  of  Equation  (29),  since 


2tt 


X 


~h 


exp[Bae  cos  E]  cosm  E ( 1  -  e^  cos^  E)  sin  E  dE  =  0 


(44 


Thus,  Equations  (30)  and  (39)  may  be  used  to  express  Equations  (23)  through 
(29)  in  the  following  compact  matrix  form 


A  e  = 


2k+n 


C  n  J n,l 


(45 
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